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Abstract 


By combining a holographic p-wave superconductor model and a holographic fer- 
romagnetism model, we study the coexistence and competition of ferromagnetism and 
p-wave superconductivity. It is found that the results depend on the self-interaction 
of magnetic moment of the complex vector field and which phase appears first. In the 
case that the ferromagnetic phase appears first, if the interaction is attractive, the 
system shows the ferromagnetism and superconductivity can coexist in low tempera- 
tures. If the interaction is repulsive, the system will only be in a pure ferromagnetic 
state. In the case that the superconducting phase appears first, the attractive inter- 
action will leads to a magnetic p-wave superconducting phase in low temperatures. 
If the interaction is repulsive, the system will be in a pure p-wave superconducting 
phase or ferromagnetic phase when the temperature is lowered. 
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1 Introduction 


In condensed matter physics, there are two kinds of critical phenomena which have at- 
tracted a lot of attention for a long time. One is ferromagnetism where the electron spins 
align to produce a net magnetization, which breaks the time reversal symmetry sponta- 
neously and happens in the ferromagnets at the Curie temperature, To (sometimes it is 
even higher than the indoor temperature). The other is superconductivity where electrons 
condense into Cooper pairs in the momentum space, which breaks the U(1) symmetry spon- 
taneously and usually happens at a much low temperature Te. Though ferromagnetism 
has been found for more than thousands years and the superconductivity was found about 
100 years ago, their physical explanations was provided after the complete quantum theory 
on materials was built. 

In a very long time, it was thought these two phenomenons are incompatible with 
each other. This is rooted in the microscopic theory of superconductivity by Bardeen, 
Cooper, and Schrieffer (BCS) [1]. In the framework of BCS theory, the superconductivity 
appears when the electrons were bounded with antiparallel spins in singlet Cooper pairs 
due to the effective attractive force coming from the lattice vibrations. When magnetic 
impurity atoms are placed in a conventional superconductor, they are capable of flipping the 
electron’s spin. Hence, impurity will suppresses the singlet Cooper pair formation, which 
causes a rapid depression of the superconducting transition temperature Te- Likewise, 
the superconductivity can screen off the magnetic field, which leads to the long-range 
magnetic order is accompanied by the expulsion of superconductivity. These features 
lead to the competition of such two orders. However, around 1980s, it was recognized 
that under special conditions superconductivity may coexist with antiferromagnetic order, 
where neighboring electron spins arrange in an antiparallel configuration. For instance, in 
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Figure 1: The experimental results on the p-wave superconducting ferromagnetic material 
URhGe. (a) Temperature dependence of the spontaneous magnetization M (T) and the in- 
verse of the magnetic susceptibility x~! (T) under the normal pressure [6]. (b)Temperature 
dependence of the resistivity of URhGe for different pressures at low temperature. The su- 
perconducting critical temperature is defined by the zero resistivity point [7]. (c) Pressure- 
temperature phase diagram of URhGe [7]. 


heavy fermion antiferromagnets, the itinerant magnetic moments have almost no de-pairing 
effect on singlet Cooper pairs, because the average exchange interaction is zero. 

The discovery of the first superconducting ferromagnet! UGeg in the year 2000 came as 
a big surprise [2]. In this material, superconductivity is realized well below the Curie tem- 
perature, without expelling the ferromagnetic order. Since then, other two superconducting 
ferromagnets have been discovered, such as URhGe [3] and UCoGe [4], which display in- 
trinsic coexistence of ferromagnetism and superconductivity. Evidence of superconducting 
ferromagnetic phase was also reported for ZrZnə in 2001 [5]. 

Figure 1 shows the experimental results of a typical superconducting ferromagnetic 
material URhGe. The ferromagnetic order is observed near To œ 9.5K, where the spon- 
taneous magnetization begins appearing and the magnetic susceptibility diverges. Below 
this temperature, the initial magnetic susceptibility is positive and the spontaneous mag- 
netization increases with decreasing the temperature. In the ferromagnetic phase, when 
temperature is decreased to near Tse ~ 0.25/ (at atmospheric pressure), the resistivity 
approaches to zero, which shows that the material is in a superconducting phase. In the 
region of T < Te, spontaneous magnetization and zero resistivity coexist, which shows 
a superconducting ferromagnetism. The P — T phase diagram shows that this kind co- 
existence can appear in a very wide pressure region (up to near 30 kbar, where the Tse 
disappears. ). 

'In this paper, we will use “superconducting ferromagnet” to denote the materials whose Curie tem- 


perature is higher than superconducting transition temperature and “ferromagnetic superconductor” to 
denote the opposite case. 


The nature of superconducting state in ferromagnetic materials is currently under de- 
bate. Early investigations [8] studied the coexistence of conventional s-wave superconduc- 
tivity with itinerant ferromagnetism. However, the scenario of spin-triplet pairing soon 
gained the upper hand [9]. For a review of phenomenological theory of ferromagnetic un- 
conventional superconductors with spin-triplet Cooper pairing of electrons, one can see 
Ref. [10]. Ref.[11] presents a general thermodynamic theory that describes phases and 
phase transitions of ferromagnetic superconductors with spin-triplet electron Cooper pair- 
ing, based on an extended Ginzburg-Landau theory. Generally speaking, the coexistence 
of ferromagnetism and superconducting is discovered in spin triplet rather than the usual 
spin singlet superconducting materials, because that p-wave pairing allows parallel spin 
orientation of the fermion Cooper pairs in unconventional superconductors [12]. This 
unconventional paired manner makes p-wave superconductivity being robust under the in- 
fluences of external magnetic field and spontaneous magnetization. So they may coexist 
with each other. 

A mean-field model for coexistence of spin-triplet pairing and ferromagnetism was de- 
veloped in [9, 13]. The model considers a uniform coexistence of ferromagnetism and 
superconductivity, i.e., the same electrons play the role for the ferromagnetism and su- 
perconductivity at the same time. Another scenario where there is an interplay between 
magnetic and superconducting orders in the same material is superconductors with spiral 
or helical magnetic order. Examples of such include ErRh4B4 and HoMogSs. In these cases, 
the superconducting and magnetic order parameters entwine each other in a spatially mod- 
ulated pattern, which allows for their mutual coexistence, although it is no longer uniform. 
Even spin-singlet pairing may coexist with ferromagnetism in this manner. 

Up to now, the theoretical investigation have been concentrated on the weak coupling 
case, where some approximations and conception of free field are still valid. The investiga- 
tions on strong correlated system in theoretical and experimental aspects have challenged 
the pictures about the materials. A crucial feature of these systems is the nonzero mag- 
netic moment of the spin-triplet Cooper pairs. However, the microscopic theory about the 
coexistence of magnetism and superconductivity in strongly interacting heavy electrons is 
either too complex or insufficiently developed to describe the complicated behavior. So it 
is still a fascinating thing to find a suitable theory to describe the coexistence and com- 
petition of the ferromagnetism and superconductivity in strong correlated system, such as 
the heavy fermion system [14], Iron-based superconductor [15] and unconventional super- 
conductor [16]. 

In the strong coupling case, the usual methods developed in condensed matter theory 
(CMT) are considered losing their efficacy. A new method, named gauge/gravity duality 
or AdS/CFT correspondence, is considered as a promising approach [17, 18, 19, 20]. This 
method relates a weak coupling gravitational theory in a (d+ 1)-dimensional asymptotically 
anti de-Sitter (AdS) space-time to a d-dimensional strong coupling conformal field theory 
(CFT) in the AdS boundary. In recent years, this duality has been extensively applied into 
condensed matter systems. Some remarkable progresses have been made in this direction. 
For example, some gravitational dual models of superfluid/superconductor [21, 22, 23], 
(non-)Fermi liquid [24, 25, 26], Josephson junctions [27, 28, 29], superconducting quantum 
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interference device [30] and magnetic properties in superconductors [31, 32, 33, 34, 35] 
have been constructed and intensively studied. The models in the AdS/CFT frame for fer- 
romagnetism/paramagnetism and anti-ferromagnetism/paramagnetism phase transitions 
have also been proposed in [36, 37]. 

In this paper, we will explore the coexistence and competition between ferromag- 
netism and p-wave superconductivity in the strong coupling case. Our tool is just the 
gauge/gravity duality. We noted that the similar topic appeared in Ref. [38], where the 
authors use two U(1) fields condense simultaneously from SU(2) model to present the 
superconducting and spontaneous magnetic orders. Our model is different from theirs. 

In order to construct a holographic model to investigate the coexistence and competition 
between ferromagnetism and superconductivity, we need a model where ferromagnetic 
properties (such as spontaneous magnetic moment, time reversal symmetry broken and 
diverged magnetic susceptibility at Curie temperature and so on) can appear independently 
from the superconductivity. This is just supplied by the model in Ref. [36]. By combining 
the complex vector field model for the holographic p-wave superconductor [39, 40, 41] and 
the real antisymmetric tensor field model for the holographic ferromagnetism in Ref. [36], 
we study the coexistence and competition of ferromagnetism and superconductivity. It 
turns out that the results depend on magnetic moment self-interaction of complex vector 
field and which phase appears first. In the case that the ferromagnetic phase appears first 
and the interaction is attractive, the p-wave superconductivity can still appear and the 
system can show ferromagnetism and superconductivity both when temperature is lower 
than a critical value. But if the interaction is repulsive, the p-wave superconductivity 
can not appear and the system will only be in a pure ferromagnetic state. In the case 
that the superconducting phase appears first, the system will show a magnetic p-wave 
superconducting phase with decreasing temperature if the interaction is attractive. If 
the interaction is repulsive, the system is in a pure p-wave superconducting phase or 
ferromagnetic phase when the temperature is lowered. 

The paper is organized as follows. In section 2, we will first describe the holographic 
model. We will give our ansatz for matter fields, the equations of motion (EoMs) and 
the expression of free energy density in section 3. In section 4, we will investigate the 
possible phases and the coexistence and competition between ferromagnetism and super- 
conductivity with different parameters. A brief summary and discussion will be given in 
section 5. 


2 The model 


In this paper, the model we are considering is the combination of the holographic p-wave 
superconductor model described by a complex vector field [39, 40, 41] and holographic 
ferromagnetism model described by a real antisymmetric tensor field [36]. The action is 


1 6 v 
S= 55 | dav [e+ te -r + Lp + Lm + Leone ’ (1) 


with 


1 V . V 
Lo = —5Plw P” — MILO + iqypu o F” — Vp, 


1 2 A 
Lu = - 704M" V „Mn — “aM Mpy — SM" Py — Vm, (2) 


Lom = —iap pM, 


where L is the AdS radius which will be set to be unity and «K? = 87G is related to the 
gravitational constant in the bulk. In the following, we will set 2x? = 1 for simplicity. A, 7 
and J are three constants with J < 0 [36], y > 0. mı, mz are the masses of the complex 
vector field p, and real tensor field M,,,. g is the determinant of the bulk metric gv and q is 
the charge of complex vector field. We define Fi, = V, Ar — VLA, and ppv = Dupo — Di py 
with the covariant derivative D, = V, — iqA,. The y’s term characterizes the magnetic 
moment of the vector field p, [39, 42]. The antisymmetric tensor M,,, is the effective 
polarization tensor of the U(1) gauge field strength Fav. Vm describes the self-interaction 
of the polarization tensor. Following Ref. [36], we take 


Vu = ZM M," M5 Mg (3) 
for simplicity. The term V, describe the self-interaction of magnetic moment of the complex 
vector field, a simple form is 5 

V, = -7 Puhr p o”. (4) 
Here O is a constant, it characterizes the feature of magnetic moment interaction of complex 
vector field. From the free energy density given later, we can see that a positive O gives an 
attractive interaction between the magnetic moment of complex vector field itself, while a 
negative value of © gives an repulsive interaction. In Refs. [39, 40, 41, 42], the term V, is 
not considered as there the aim is to study the superconductivity of the model. However, 
we can see that this term will play an important role in the following when we consider 
the magnetic properties of the model. 
One may find that the Lagrangian £m has a little difference from the form in [36]. This 
is just for convenience. In order to simplify our discussion in the probe limit, we make 
following transformations 


My > AM, J > A?J, pu > àp a Q/A, O 3 OA. (5) 
Under these transformations, action (1) can be rewritten as 
6 
S= par r + TaT Fp F” +X (L, + Lum + Lom) l, (6) 


with 
1 HAJUT ms HV 1 Hv 
Lm = =N M V Myr _ q” Mw m a Puy — Vu (7) 
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and the others are kept the same as the forms in (2). The probe limit corresponds to 
keeping all the quantities finite with A —> 0. In this limit, we can fix the background 
geometry and Maxwell field and only consider the dynamics of complex vector field and 
polarization field. 


3 EoMs and free energy density 


3.1 Ansatz and EoMs 


As in Ref. [36], we will work in the probe limit with A — 0, by which we can fix the 
geometry background and neglect the back reaction of matter fields on the background 
geometric and Maxwell field. The equations for complex vector field and polarization field 
read 


D" pup — Mi py + tgp” Fi, + Oppo — tap’ M,, = 0. 


V?2M,, —m2M,, — JM,,°M;"M._.,, — 2i Lf. = (8) 
pw 2th aw pts Mery AP (uP) ww = O. 


We take the AdS Reissner-Nordstr6m (RN) black hole with a planar horizon as the back- 
ground metric [43] 


2 2 2 2 2 d? 
(9) 


2 
f(ir)=1- 3 P A, = p(1 — 1/r)dt. 


Here the horizon radius has been set to r, = 1 and u can be identified with the chemical 
potential in the dual field theory. The temperature of the boundary theory is 


o1 
At 


T (3 — u’). (10) 


A self-consistent ansatz for complex vector field and polarization field is 
Mw = —p(r)dt ^ dr + h(r)dx Ady, py = prdx + ie pydy. (11) 


As in Ref. [39], we can take p,(r) to be real by some suitable U(1) gauge. However, the 
y-component of p, then can not be set to be real. So we have to assume that p, and 0 
are real functions depending on r. Taking the ansatz into equations (8), we find that only 
when e(") = +1 for Vr € (rn, 00) can we find a self-consistent solution. Without loss 
generality, we can assume (r) = 0. Considering the symmetry between p, and py, we 
assume that 


Py = C(") Px (12) 


in what follows. Thus we can reach the following equations for the components of matter 
fields 


; : Jh 2f 4 2ca p? 
Lo arek. Pp 
P 2 P Oe m? cha 
" (eae / j = o ee  — Y, 13 
Pat (ot dat aaa Gag pe ree (13) 
/ 2 WN, -A(O tah 
+ (F424 7B) ae ao 


The equation for p(r) decouples from the above equations, so we will not write down it 
here. 
Near the AdS boundary, the linearized equations give following asymptotic solutions ? 


p®-1)/2 -(61+1)/2 


ô 
F prf C= GTI +C, 


h(r) = h, r@+®2)/2 a h_r-9/2 


where 6; = \/14+ 4m? and 62 = \/17+4m3. According to the AdS/CFT dictionary, 
Pz+, Cy and h, are sources terms for corresponding operators, while p,_, c_ and h_ are 
vacuum expectation values, respectively. As in Ref. [36], we need impose the condition 
h} = 0 for the polarization field and pr} = c+ = 0 for the complex vector field. This is 
consistent with the spirit of AdS/CFT correspondence: one requires that the condensation 
and magnetization happen spontaneously. 

At the horizon, we impose regular conditions for p,, c(r) and h, which give following 
relationships for the initial values at r = rp 


Px = Pat (14) 


d 


1 
Pr = =a ho(mj + cah — cop), 


ATT 
aye 2 
nf — op — POR — (15) 
1 
! = — (1 -e ? + ah). 
č nT! e) (cO + ah) 


Note that there exists a symmetry as {ps — Py, Py — Px}, by which we can set that 
lPy(7n)| < [Px (rn) |, 1.€., 
—1<e(ra) <1. (16) 


Thus once given the value of parameters {9, a, J, q, m?, m3, T} and some suitable 


initial values {h(rp), Px(Tn), C(rn)}, we can integrate equations (13) to obtain the whole 

solutions matching the boundary conditions hy = pz, = c+ = 0 at the AdS boundary. 
Our numerical results show that °, except for case that c(r,) = h(r,) = 0 or e(r) = £1, 

the integration will meet a divergency at somewhere in r, < r < oo and the equations (13) 


?The asymptotic solution of c(r) depends on the source free condition of py. When pr+ 4 0, asymptotic 
solution of e(r) is c = c} +c- rÀ. 

3We use the function ode45 in MATLAB R2012b with the relative and absolute errors 10713 to solve 
equations (13) numerically. 


do not have physical solutions. When c(r;,) = h(r;,) = 0, according to the relationship 
of initial values in (15), we can find that e(r) = h(r) = 0, which is just the p-wave 
superconductor solution found in Ref. [39]. When c(r;,) = £1, the solution for c is e(r) = 
+1, i.e., Pg = py, which leads that the nontrivial solutions are either pure ferromagnetic 
phase or p-wave superconductivity phase with nonzero magnetization if a # 0. From the 
equations (13), we can see that the transformation a — —a is equivalent to fix a but make 
the transformation of c(r) => —c(r). So the case of a < 0 is equivalent to the case of 
a > 0 with exchanging the results of c(r) = +1. In the special case that a = 0, from the 
equations (13), we can see that the equations for h and py are decoupled with each other. 
The ferromagnetic phase and p-wave superconducting phase have no interaction with each 
other +. We are here not interested in that case. In the following sections, therefore we 
will set a > 0. 


3.2 Free energy and magnetic moment 


Let us now compute the on-shell action and give the thermodynamic of the dual boundary 
theory. Here, we use grand canonical ensemble by fixing the boundary chemical potential. 
For convenience, we can set u = 1 in the numerics. The results for other chemical potential 
values can be obtained by scaling relations. We will first fix the horizon radius r, to 
solve the equations (13) and compute the free energy density, then we use the scaling 
transformation to obtain the results in the grand canonical ensemble. In gauge/gravity 
duality, the Gibbs free energy F can be obtained by temperature timing the on-shell bulk 
action with Euclidean signature. Since we work in the probe approximation, we can ignore 
the gravity part. Given that the system is stationary, the Euclidean action is related to the 
Minkowskian one by a total minus. Using the equations of motion (8) and the source free 
conditions, we can get the free energy density contributed by vector field and polarization 
field as = 

2 = x | dr./—g (-¥, + Luwr, = Vm ~~ aa baT) , (17) 

V ee 4 2 
where V is the area spanned by coordinates x and y on the boundary. Taking the 
ansatz (11) into the free energy density, it turns out to be 


F oo ht 2 d 2 
y=% f ar (25 - A t), (18) 


4r6 r2 r2 


Note that the contribution to the free energy density from p(r) is not relevant to our 
discussions, therefore here we have neglected that part in (18). The integration is finite 
when pr} = h} = 0 and m? > —1/4, m3 > —4. In addition, note that there are two 


symmetries in the free energy and EoMs (13) such as 


{Pz > —Pz, Py > —py}, {c > —c, h > —h} (19) 


4This is the consequence of the probe limit. Once the back reaction is taken into account, even when 
œ = 0, p, and Mpy will interact each other through gravity background. 
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which make we can specify p,(r,) > 0 and h(r;,) < 0. It is a general requirement that the 
functions of p, and h do not have zero points in the region of |r}, o0). Thus we can assume 
px(r) > 0 and A(r) < 0 in the region [rp,00). According to the definition of magnetic 
moment density and using the expression (17), we have 


1 OF = i 
N = —— li — = — — d 20 
V B20 Ge [ D2” ay) 
which is the same as in Ref. [36]. By the dictionary of AdS/CFT, the expectation value 
of p-wave superconducting order parameter is a complex vector P, whose mode is P = 


V1+4+c?|p,_|. Here it is worthwhile to mention that though the complex vector field does 
not appear in the final expression of magnetic moment density, it makes contribution 
through the mixture terms in equations (13): once py and c do not vanish, h will not 
vanish. 

In the pure p-wave model, when nontrivial solutions of p, appear, the global U(1) and 
spatial rotation symmetries are broken spontaneously. In this model, it is also true. In 
addition, there is an another possible symmetry breaking in this model. If one notes the 
fact that external magnetic field B will be transformed in —B under the time reversal 
transformation, by the expression of free energy density in (17), we have following rules 
for time reversal transformation, 


h> =h, Py > =Py (21) 


So when A # 0 or py = +p, # 0 (they both lead to nonzero magnetic moment), the time 
reversal symmetry is broken spontaneously. This agrees with the fact that a spontaneously 
magnetized phase will break time reversal symmetry spontaneously. 

From Refs. [39, 36], we can see that the complex vector field and polarization field can 
both condense in low temperatures in an AdS RN black hole background. We take T;,9 
and Toco as the critical temperatures of p, and h, when a = 0, in the AdS RN black hole 
background, respectively. With decreasing the temperature of the system, three interesting 
questions appear immediately: 

(1) If Too > Tsco, then the system will enter into the ferromagnetic phase first. Can 
a p-wave superconducting phase still appear at lower temperature? If yes, is the critical 
temperature still T’,.9? 

(2) If Tso > Too, then the system will enter into the p-wave superconducting phase 
first. Can a ferromagnetic phase appear at some lower temperature? If yes, is the critical 
temperature still Teo? 

(3) In an enough low temperature, for example, near the zero temperature limit, can 
the p-wave superconductivity and ferromagnetism coexist in this model? 

In what follows, we will consider all these questions. We will first solve equations (13) 
with the regular condition (15) at the horizon and source free conditions at the AdS 
boundary. In general, depending on parameters, the system has four kinds of solutions: 
one is a trivial solution without vector and tensor hairs, one is a solution with vector hair, 
but no tensor hair (which describes the pure p-wave superconductivity phase), one is the 
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solution with tensor hair, but no vector hair (which describes the pure ferromagnetic phase), 
and final one is the solution with both hairs (which describes both the superconductivity 
and ferromagnetism coexistence phase). The phase is physical favored if it has the lowest 
free energy. 


4 Coexistence and competition 


4.1 Superconducting ferromagnet 


Let us first consider the case with Teco > T'sco, i.e., the ferromagnetic phase appears first. 
According to the equation for c in equations (13), we find c 4 0 when h Æ 0. So there is not 
a phase such that {h < 0, pz #0, py = 0}. When we decrease the temperature to be lower 
than Tco, five kinds of solutions may appear. They are phase A {h = ps = py = 0}, phase B 
{h <0, pe = py = 0}, phase C {h = py = 0, pa Z 0}, phase Dı {h < 0, pz = py + 0} and 
phase Də {h < 0, ps = —py # 0}, which corresponds to normal phase, pure ferromagnetic 
phase, pure p-wave superconducting phase and two kinds of superconducting ferromagnetic 
phases, respectively. 
For small p,, we can get its effective mass square from equations (13) as 


(22) 


In general, h(r) does not have zero point in the region of r, < r < oo in the condensed 
phase. With the choice h < 0, a > 0, we have ah < 0. If c = 1, the effective mass square of 
complex vector field will decrease by the condensate of h. If c = —1, instead the effective 
mass square of complex vector field will increase by the condensate of h, which leads the 
critical temperature of complex vector field to decrease. In this case, if —ha is enough 
large, the complex vector field can not condense even in zero temperature. These imply 
that when we decrease the temperature, the instability of complex vector field will appear 
in the manner of py = pz Or Py = —Pz for small a and in the manner of p, = p, for enough 
large a. This instability tells us that in the ferromagnetic phase, the marginally stable 
mode for complex vector field can still appear at some temperature less than Teo. The 
analysis can also be confirmed directly by solving equations (13) numerically under the 
case of neglecting the p2 terms, i.e., 


í dhe 27? 4 m 
h! + =R a ee E E |) 
ikl ae aa) a 


23) 
/ 242 2 ( 
n" f 2, , go mi cha 

a — Ht ) oe =0, 
aaa f sae a E fr? ot) 
with the initial value of p,(r,) = 1. As a typical example, we choose parameters as 
m? = —3/16,m3 = —3, J = —1 and q = 1.3 (for other parameter values, the results are 


qualitatively similar) . The results are shown in figure 2. In the left plot, we show p4 as a 


11 


0.05 


—c=l 


—c=-1 


-0.05} = . ; Too 
‘Teo: | 
Ta i : 
2 Tex 
-0.1 nE ei a 6 , ! a _ C To 
0 0.2 0.4 0.6 0.8 1 12 0 02 04 0.6 08 1 
T / Too a 


Figure 2: Left: p, as a function of temperature T when h < 0. The green line and blue 
line stand for c = —1 and c = 1, respectively. The red dashed line stands for the case in 
the pure AdS RN background with h = 0. Right: Tsc/Tco and Tse /Tco as functions of a. 
We take m? = —3/16,m?2 = —3, J = —1,q = 1.3. Here Teo œ 0.00925, Teo = 0.8135Tco 


function of temperature T for c = +1 in the case of a = 0.1 when h < 0. The green line 
and blue line stand for c = —1 and c = 1, respectively. The red dashed line stands for the 
case in the pure AdS RN background with h = 0. In the case of c = 1, there is a critical 
temperature Tse less than Teco but higher than Tsco to make p4 = 0. In the case of c = —1, 
there is a critical temperature Tsw less than Tsco to make p} = 0. In the right plot, we 
show T,e/Tco and Tse /Tco as functions of a. When we increase a, Tse will increase but 
T,» decreases. Numerical results show that there is a critical value at a = a, ~ 0.57649, 
less than which Tse will be less than zero. In appendix A, we give the method to compute 
the critical value a. 

Note that equations (23) and figure 2 only show that there is a marginally stable mode 
for the complex vector field at T = T,,. or T = 7... When temperature is below these 
critical temperatures, whether the complex vector field can condense and which one of two 
phases {h < 0, ps = —py # 0} and {h < 0, px = py # 0} can appear in the physical phase 
space are determined by their free energy density. In order to find the phase diagram, we 
have to solve equations (13) numerically to compute the free energy of possible solutions. 
It turns out that the results depend on the sign of ©. The possible phases and the physical 
favored phase in different temperature regions are summarized in table 1. 

In the case of O > 0, we find that the the phase D; is physical favored when T < Tze- 
This means that the system displays ferromagnetism and superconductivity both in low 
temperatures. In addition, by computing the ferromagnetic order parameters and p-wave 
superconducting order parameter, we see they are both continuous at two critical temper- 
atures. With lowering the temperature, the system will first transit into the ferromagnetic 
phase at T¢o and then into the superconducting ferromagnetic phase at Tse through two 
second order phase transitions. The critical temperature of complex vector field grows 
when a gets increased. This shows that, in the ferromagnetic state, the interaction be- 
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Phases in the case of a < a, 


Te STE 


è A, D: A, B, dD, 
ou? < Ta) | bud 


PhialOso| A |B | 
Physical <0)[ A [SSS 


Phases in the case of a > a, 


T 
oee A [AB [ABD CET eTa) 


Physicel@>0)[ A |B | So 
Physical <0) APB 


Table 1: The possible and physical phases in the case of Teco > Tso. Phase A is {h = py = 
Py =O}. Phase B is {h < 0, pz = py = 0}. Phase C is {h = p} = 0, pz #0}. Phase D, is 
{h <0, pr = py #0}. Phase Dg is {h <0, pz = —py F 0}. 


tween p-wave pairs and spontaneous magnetic moment will promote the appearance of 
p-wave superconductivity. On the other hand, in the case of © < 0, though the solutions 
of pr # 0 exist, they are not physical favored because they have higher free energy than 
the solution of py = 0. As a result, in this case, there is only a ferromagnetic phase (phase 
B) when T < Too. 

Here we only show the example for the case of a < ae. The case of a > a, is similar 
except for the difference that the solution for the phase Də does not occur. In the left plot 
of figure 3, we show the free energy density of phases Dı and D» in the case of O = +1 and 
phase B with respect to temperature. Note that here the on-shell free energy density of 
phases A and Č is zero in the probe approximation. One can see that phase D; and phase 
B are physical favored when © = +1 respectively when T < Te- Thus we see that the 
physical favored phase depends on the sign of O. In the case of O > 0, with decreasing the 
temperature, the system will first goes into the ferromagnetic phase when T < Teco and then 
into the p-wave superconducting ferromagnetic phase when T < Tse. However, there are 
two superconducting ferromagnetic phases D,({p, = py}) and Dj({p, = —p,}). Which has 
a lower free energy than the pure ferromagnetic phase when T < Te? It turns out that the 
physical favored one is phase Dı. The numerical results show that the ferromagnetism and 
superconductivity can coexist in the whole region of T < Tse. It seemingly indicates that 
the ferromagnetism and superconductivity can coexist even in the zero temperature limit. 
However, in the case of © < 0, we see that although the superconducting ferromagnetic 
phases exist (phase Dı and D2), these two phases have higher free energy density than the 
pure ferromagnetic phase. In this case, therefore the physical favored phase is the pure 
ferromagnetic phase and the p-wave superconductivity will not occur. 

In the right plot of figure 3, we show the superconducting order parameter P and 
spontaneous magnetic moment density N as functions of temperature T in the case of 
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Figure 3: Left: The free energy density with respect to temperature in the case of © = 
+1. Right: The p-wave superconducting order parameter P and spontaneous magnetic 
moment deneny N as functions of a in the case of © = +1 in the physical favored 
phase. Here m? = —3/16,m32 = —3, J = —1,q = 1.3 and a = 0.1. A_=1+(1+6,)/2. 
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Figure 4: The p-wave superconducting order parameter P and spontaneous magnetic mo- 
ment density N with respect to temperature near the critical temperature in the Phase D4. 
Numerical fittings ewe that Pu-4- œ 0.3721,/1 — T/T. and N/p œ~ 0.1588,/1 — T/Too. 
Here m? = —3/16,m3 = —3, J = —1,q = 1.3 and a = 0.1. AL =1+(146,)/2. 
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© = +1 in the physical favored phase. We see that there is a second order phase transition 
at T = Too, below which the system enters into a ferromagnetic phase. When temperature 
decreases below Tse, the situations depend on the sign of ©. In the case of O > 0, there is a 
second order phase transitions at T = Tse such that the system will transit into the p-wave 
superconducting ferromagnetic phase. From the curve of N in the right plot of figure 3, 
one can see that the condensation of p-wave order will increase the magnetic moment. This 
implies that the p-wave pair carries a nonzero magnetic moment, which contributes to the 
total magnetic moment of the system. In figure 4 we plot the behaviors of the magnetic 
moment and p-wave order parameter with respect to temperature in the phase Dı near 
the critical temperature, which clearly shows a square root behavior for both quantities. 

Let us make a brief summary for this subsection. In the case of Too > Tsco, i.e., the 
case with the ferromagnetic phase appearing first with decreasing the temperature, whether 
the p-wave superconductivity can appear depends on the sign of ©. If © > 0, there is a 
critical temperature Tse which is lower than Too but higher than T,.9. When T < Tse, the 
p-wave superconductivity can appear and the system will show the ferromagnetism and 
superconductivity both. Even in the near zero temperature limit, they can coexist. On 
the other hand, if © < 0, the p-wave superconductivity can not appear and the system 
will only be in a pure ferromagnetic phase. These results are summarized in table 1. 


4.2 Ferromagnetic superconductor 


Let us now consider the other case that 7.9 > Too, i.e., the case where the p-wave super- 
conducting phase appears first. When Too < T < Tæ, according to the equations (13), 
there may exist three kinds of p-wave superconducting phase. One is just the usual p-wave 
superconducting phase C({h = py = 0, ps # 0}), the other two are new superconducting 
phase Fy with {A < 0, pr = py # 0} and Ey with {h < 0, pr = —p, # 0}. Thought phase 
E, and E> have nonzero magnetic moment, they are different with phases Dı and Dap, 
because the appearance of nonzero h in the former two is induced by the p-wave pair but 
in the latter two is spontaneously produced. Thus we have two questions as follows. Can 
these three solutions exist ? And which one is physical favored when temperature is in the 
region Too < T < Tyco? 

Our numerical results show that the answers depend on the sign of ©. If O > 0, 
the phases A, C and EF, can exist and the phase EF, is physical favored which has the 
lowest free energy. Therefore the system will show magnetism once it goes into the p-wave 
superconducting phase. This case is very similar to the Anderson-Brinkman-Morel (ABM) 
phase in *He superfluid [44], where the superfluid phase is also of magnetism. Though 
the magnetism and superconductivity appear together, it has an essential difference from 
phases D; and D» just as we mentioned before. So the phase Æ; (and EF) should be called 
“magnetic superconducting” phase rather than ferromagnetic superconducting phase. This 
difference can also be shown in the magnetic moment density near the phase transition 
point. In the phases Dı and D2, we see that the N shows a square root behavior with 
respect to 1—T/Too (see figure 4), while in phase Æ, we see that N has a linear relationship 
with respect to 1 — T/T (see figure 6). On the other hand, if © < 0, we find that the 
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Figure 5: Left: The free energy density with respect to temperature in phase F and 
phase B. Right: The p-wave superconducting order parameter P and induced magnetic 
moment density N with respect to temperature in phase E,. Here m? = —3/16,m3 = 
—3, J = —1,q = 1.4,0 = 1 and a = 0.1. A- = 1 + (1 + ô1)/2. 


phase Æ; and E» do not appear and the equations (13) only have the trivial solution and 
the pure p-wave superconductivity solution. As a result, in this case, the system is in a 
pure p-wave superconducting phase without magnetism in the region of Too < T < Tico. 
As an example, let us consider the parameters as m? = —3/16,m3 = —3, J = —1, q = 
1.4 (again, the results are qualitatively similar for other parameter values, the only re- 
quirement is to have T..9 > Too). In this case, we have T,.9 ~ 1.83837 ¢9. The free energy 
density, the p-wave superconducting order parameter P and induced magnetic moment 
density N with respect to temperature in the phase Æ are shown in figure 5. Because the 
on-shell free energy for phases A and C is zero in the probe approximation, we see that 
phase E; has the lowest free energy. Furthermore from figure 6, we see that the p-wave 
superconducting order parameter P have a square root behavior with respect to 1 —T/Tsco 
but the induced magnetic moment density has a linear relationship with 1 — T/T ..9 when 
TST g 
When temperature is lower than Too and © > 0, phases B can also appear. The 
numerical results show that the phase F is still the lowest free energy phase (see figure 5 
as an example). So there is no phase transition at T = Too. If O < 0, we have known that 
when T > Too, this system is in a pure p-wave superconducting phase C with {h = p, = 
0, pz # 0}. At T = Too, because the equation for h is a homogeneous one (py = CPx, in 
phase C we have c = 0), there is critical point for h. However, the equation for c(r) in (13) 
restricts c(r) to be 0 or +1. When h ¥ 0, the only solution for c(r) is e(r) = +1. This 
indicates that in this case the solutions for (13) are either just phase FE, (or E2) or phase 
B, i.e., phases {h < 0, ps = tp, # O} or {h < 0, ps = py = 0}. However, our numerical 
calculations show that such solutions {h < 0, ps = +p, # 0} do not exist. Thus the possible 
phases are only phase A, phase B and phase C. In this model, the on shell free energy is 
zero for phase A and C but is negative for phase B. So the phase {h < 0, ps = py = 0} 
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Figure 6: The behaviors of N and P and near the critical temperature To in the phase Æ. 
Numerical fittings show that Pu~4- ~ 0.3699,/1 — T/T, and N/u ~ 0.0294(1 — T/T). 
Here m? = —3/16,m3 = —3, J = —1,q = 1.4,0 = 1 and a = 0.1. A- = 1 + (1 + ô1)/2. 


Phases in the case of O > 0 


C Possible | A | AEC [A ECB 
aa a a 


Phases in the case of O < 0 


[Physics [A [oc _| B_ 


Table 2: The possible and physical phases in the case of Teo > Too. Phase A is {h = pz = 
Py = 0}. Phase B is {h <0, ps = py = 0}. Phase C is {h = py = 0, ps #0}. Phase F is 
{h <0, pc = py #0}. Phases Hy is {h <0, pz = —py 4 0}. 


is physical favored in the case of © < 0 when T < Teco. The ferromagnetism can still 
appear from the p-wave superconducting phase at the same critical temperature To, but 
the p-wave superconducting phase will disappear. In other words, the superconductivity 
and ferromagnetism can not coexist in the case of © < 0. In figure 7, we plot the difference 
of free energy between phase C and phase B and magnetic moment density N with respect 
to temperature in phase B, where F = Fyhase B — Fphase c, from which we can see that 
the pure ferromagnetic phase has lower free energy than the pure p-wave superconducting 
phase. 

Let us make a brief summary for this subsection. In this case of Teo < Tsco, i.e., the case 
that the superconducting phase will appear first, the results depend on the sign of ©. If 
© > 0, the ferromagnetic phase can not appear but the magnetic p-wave superconducting 
phase will appear in region of T < To. If O < 0, the system will be in a pure p-wave 
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Figure 7: Left: The difference of free energy between phase C and phase B. Here ôF = 
Fphase B — Fphase c. Right: The magnetic moment density N with respect to temperature 
in phase B. Here m? = —3/16,m3 = —3, J = —1,q = 1.4,0 = —1 and a = 0.1. 


superconducting phase in the region Teo < T < Tsco and a pure ferromagnetic phase when 
T < Too. All the results are summarized in table 2. 


5 Summary and discussions 


In this paper, by combining the complex vector field model for the holographic p-wave 
superconductor and the real antisymmetric tensor field model for the holographic ferro- 
magnetism, we have investigated the coexistence and competition of ferromagnetism and 
superconductivity in the holographic setup. Depending on model parameters, we found 
that the model shows rich phases in low temperatures. The study is done in the probe 
limit, the background geometry is taken to be an AdS RN black hole with a planar horizon. 

In the case of Tog > Tyco, i.e., the case where the ferromagnetic phase appears first, 
whether the p-wave superconductivity can appear depends on the sign of ©, the interac- 
tion strength of magnetic moment of the complex vector field. If © > 0, there is a critical 
temperature Tse which is lower than Teco but higher than 7-9. When temperature is higher 
than Te, the system only shows the ferromagnetism. When T < Tse, the p-wave super- 
conductivity can appear and the system will show ferromagnetism and superconductivity 
both. Because of the the spontaneous magnetization, the critical temperature of p-wave 
condensation is higher than the critical temperature without the ferromagnetic phase, and 
increases with the increasing of interaction strength between complex vector field and an- 
tisymmetric tensor field. Even in the near zero temperature limit, the magnetism and 
superconductivity can coexist. But if © < 0, the p-wave superconducting state can still 
exist but it is not the lowest free energy state. So the superconductivity can not appear 
and the system will only be in a pure ferromagnetic state. 

In the case of Too < Teco, i-e., the case where the superconducting phase appears first, 
the results also depend on the sign of ©. If © > 0, in the region of Tog < T < Teo, the 
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system will show the p-wave superconductivity and a kind of induced magnetism. The 
superconductivity and magnetism appear both, however, it is a magnetic superconducting 
phase rather than a ferromagnetic superconducting phase, because the magnetic moment 
is not spontaneously produced. The magnetic moment is proportional to Teo — T rather 
than /7-o0 — T near the critical temperature. When temperature is lower than Teco, the 
ferromagnetic phase B can exist, but it has higher free energy than phase Æ. So in the 
whole region of T < Tsco, the physical favored phase is magnetic p-wave superconducting 
phase Æ. On the other hand, if © < 0, when temperature is less than Teo, the system 
will be in the pure p-wave superconducting phase without magnetism. If temperature is 
lower than Too, the system will transit into the pure ferromagnetic phase from the pure 
p-wave superconducting phase. Therefore the ferromagnetism and superconductivity can 
not coexist in the case of O < 0. 

Now let us discuss some implications of our results. We have seen that the sign of 
© plays a crucial role in this model. This phenomenological parameter in (4) describes 
the self-interaction between the magnetic moments of complex vector field. A positive O 
means that the magnetic moments with same direction are attractive, while a negative 
O indicates that the magnetic moments with same direction are repulsive. From tables 1 
and 2, we can see that the ground state in the near zero temperature limit only depends 
on the sign of O. If we translate these attraction and repulsion into the boundary theory, 
then our results can be understood well. Since p-wave pair is spin triplet, it can be in the 
state of spin-up or spin-down. So every p-wave pair carries magnetic moment of +2up. 
If © > 0, which means that the p-wave pair will attract the pair which has the same 
magnetic moment direction and repulse the one which has opposite magnetic moment. 
So under the influence of spontaneous magnetization, the p-wave pair will be enhanced 
and survive. In addition, the magnetic moment of p-wave pair will tend to align along 
the direction of spontaneous magnetization, which increases the total magnetic moment 
of the system. As a result we indeed see the ground state is the phase where the p-wave 
superconductivity and ferromagnetism coexist. However, if © < 0, the p-wave pair will 
repulse the pair which has the same magnetic moment direction. So in the region where 
superconductivity dominates, the p-wave pair will align without net magnetism and the 
system is in a pure p-wave superconducting phase. When T < Too, the ferromagnetism 
will appear. Under the influence of spontaneous magnetization, the magnetic moment of 
p-wave pair will be compelled to align the same or opposite (depends on the value of a) 
direction of spontaneous magnetization, which leads to the magnetic moment of p-wave 
pair has same direction. But the p-wave pairs which have same magnetic moment directions 
will repulse each other, so the p-wave pair is not stable and will be de-paired. Thus the 
system can only be in the ferromagnetic phase. 

Finally let us make some additional comments on superconducting ferromagnetic ma- 
terials, since the Curie temperature is higher than superconducting critical temperature 
in general. Although superconductivity in ferromagnets was predicted more than 30 years 
ago, it took many years before the first material UGe2 was discovered and the research 
in superconducting ferromagnets has just begun recently. The main reasons why we say 
superconducting state in superconducting ferromagnets is unconventional are (i) Cooper 
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pairing carry magnetism and (ii) the gap structure of superconducting has a lower symme- 
try than the crystal lattices [45]. Let G represent the point-group symmetry of the lattice, 
T denote time reversal symmetry, and U(1) be the gauge symmetry. In the paramagnetic 
state (T > Too > Tsc) the symmetry group is given by G x T x U(1). In the ferromagnetic 
phase (T < Too) time-reversal symmetry is broken, and in the superconducting phase 
(T < T,. < Teo) gauge symmetry is broken as well. The coexistence of such two critical 
phenomena offers an attractive playground for the investigation of new phenomena, like 
the elusive spontaneous vortex lattice, the influence of spin-triplet superconductivity on 
the ferromagnetic domain size, control of tunneling currents by magnetization and so on. 
Also, it is a central issue in the understanding of superconductivity itself by the interplay of 
magnetism and superconductivity. Research on superconducting ferromagnetic materials 
will help us to expound how magnetic fluctuations can arouse superconductivity. This fun- 
damental insight might turn out to be crucial in designing new superconducting materials 
with high transition temperatures. We hope that the holographic model in this paper or 
the correspondence of AdS/CMT can give some helpful guidance in future. 
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A The method to compute a, 


In this appendix, we will give the method to compute the value of œe. We need to solve 
following equations in the zero temperature case 


/ 2 / 2 
we Eats (2-2 -Nao 
r 


(24) 
" F 2 ; 2p? 2 h 
Here 
f(r) =1-4/r2 +3/r*, d(r) = V3(1 —1/r). (25) 


By analyzing the behavior of h and p,(r) near the horizon, we find the regular form for 
them are, 


h(r) = ho+h(r= 1) +, pe=(r —1(1+ pa(r—1)+--*). (26) 
Substituting them into equations (24), the leading orders of the solutions give the coeff- 


cients in (26) as 
/ 1 1 


1 1 Be 
= —~+-—,/9—6ah 6m? — 3q?, pri = — : 
Bo a+ gV aho + Om; qf, Px Bo +1 
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(27) 


Taking (27) and (26) , we can integrate equations (24) from the horizon to the AdS bound- 
ary and get the solutions outside the horizon. When the values of m7, m3, q and J are 
given, we can treat hı and a as shooting parameters to match the boundary conditions 
h} = Pet = 0. Since the equation for h is dependent of pz, we can solve h first. Numerical 
result shows that hy ~ 15.40714197. Then we put the solution h into the equation of pz 
and treat a as a shooting parameter to match the boundary condition p+ = 0. Finally 
we find that there is solution with a, ~ 0.57649. 
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